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Due to the phase interference of electromagnetic wave, one can recover the total image of one 
object from a small piece of holograph, which records the interference pattern of two laser light 
reflected from it. Similarly, the quantum superposition principle allows us to derive the global 
phase diagram of quantum spin models by investigating a proper local measurement. In the present 
paper, we study the two-site entanglement in the antifferomagnetic spin models with both spin-1/2 
and 1. We show that its behaviors reveal some important information on the global properties and 
the quantum phase transition of these systems. 
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I. INTRODUCTION 

Quantum entanglement in the ground-state of stro ngly 
correlated systems [ HBI5 | S | ^Ib | ^ |SI5|^| iiIIi3 . KL3 . 

El 13 13 111 13 13 liolillillia lillialijhas been 

intensively studied in recent years. Its non-trivial behav- 
ior in these system at quantum phase transition point 
[27l | attracted many physicists' interest. Most of previ- 
ous works focused on the spin models. And the concur- 
rence, a measure of entanglement of two qubits |28| , has 
been widely used in studying these systems. For exam- 
ple, Osterloh et.al studied the concurrence between 
two spins located on a pair of nearest-neighbor sites in 
the transverse-field Ising model Q . They found that this 
quantity shows singularity and obeys the scaling law in 
the vicinity of the quantum phase transition point of the 
system. On the other hand, for the XXZ model, the 
concurrence is a continuous function of the anisotropic 
parameter and reaches its maximum at the critical point 

However, the concurrence is very short-ranged. It van- 
ishes quickly as the distance between two sites increases. 
For a substitution, Verstraete et al 0] proposed the con- 
cept of localizable entanglement. They showed that it is 
long-ranged in quantum spin systems |18l Il9l |20j . Al- 
ternatively, studies on the block-block entanglement be- 
tween two parts of the system have established a connec- 
tion between conformal field theory and the critical phe- 
nomena in the condensed matter physics pH l22j . Fur- 
thermore, since real systems consist of itinerant electrons, 
some authors generalized this concept to the entangle- 
ment in lattice fermion systems 0, HH Hfl • For exam- 
ple, for the extended Hubbard model, we showed that its 
global phase diagram can be identified by the local entan- 
glement [24j. Therefore, one is convinced that the entan- 
glement of the ground state, like the conductivity in the 
Mott-insulator transition and quantum Hall effect, and 
magnetization in the external-field-induced phase transi- 
tion, plays also a crucial role in quantum phase transi- 
tions. 

In classical optics, by recording interference pattern of 



two laser reflected from one object, its global informa- 
tion can be recovered from a piece of holograph due to 
the classical interference. Similarly, the reduced density 
matrix of a quantum system contains not only informa- 
tion on its subset, but also the correlation between this 
subset and the rest part of the system. This fact al- 
lows us to study some global properties of the system by 
investigating a part of it, as we can see from the non- 
trivial behavior of block-block entanglement [2ll [2^ in 
the quantum spin models, as well as the local entangle- 
ment in extended Hubbard model [24[. However, for the 
block-block entanglement, the size of the reduced density 
matrix, which is required for calculation, grows exponen- 
tially as diameter of the block increases. Obviously, it 
is impossible and, also, unnecessary to consider the case 
in which the size of block is comparable to macroscopic 
length. On the contrary, we would like to derive some 
properties of the system in a simple way. Therefore, we 
shall consider the block which contains sufficient infor- 
mation to reveal some global information of the system. 
Local entanglement, as a limiting case of block-block en- 
tanglement, plays such a role. 

The main purpose of this paper is to study the entan- 
glement between a local part and the rest part of the 
system in the ground-state of quantum spin models for 
both spin S — 1/2 and 1. Technically, to show the change 
of symmetry in the ground state of the system, we need 
to investigate the entanglement between a pair of neigh- 
boring spins and the rest part of the system. It is due to 
the fact that the ground states of these models are usu- 
ally spin singlet state. Our intention is to show that the 
global properties of the system can be also well under- 
stood from the behavior of local part due to the quantum 
coherence. Although the local part consists only of two 
spins, we find that it is sufficient to describe the quan- 
tum phase transition in these systems. Furthermore, we 
show that the change of symmetry in the ground state 
at the isotropic antiferromagnetic transition point yields 
a maximum value of local entanglement. On the other 
hand, at the ferromagnetic transition point, the singular 
behavior of the local entanglement can be clarified from 



the point view of infinite degeneracy and level-crossing. 
In the two-dimensional case, the cusp-like behavior of 
the local entanglement at the transition points implies 
the existence of long-range correlation in the thermody- 
namic system, which is absent in the one-dimensional 
case. 



II. ONE-DIMENSIONAL SPIN-1/2 XXZ MODEL 

First, we consider the spin-1/2 system. The cor- 
responding Hamiltonian of the antifcrromagnetic XXZ 
model reads 
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where S*f, S? and 5f are spin-1/2 operators at site i 
and A = J z /J x (J x = J y ) is a dimensionless parameter 
characterizing the anisotropy of the model. The sum is 
over all pairs of nearest-neighbor sites i and j. We impose 
the periodic boundary condition on the system. Thus, 
the choice the neighboring spin pairs is independent of 
site index. 

It is not difficult to prove that the Hamiltonian com- 
mutes with the z-component of total spins S^ otal = 
Y^i^i- Therefore, the eigenstate of the Hamiltonian is 
also the eigenstate of the z-component of total spin op- 
erator S^ otal . As a result, the reduced density matrix of 
the neighboring spins has the following form 
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with respect to the standard basis | TT), | }[), \ []), \ [[). 
Furthermore, the elements in the density matrix p/m can 
be expressed in terms of the correlation functions. We 
have 



1 



z = {S?Sf) + (S?S?), 
Wl = W2 = \- (S?S[ 



(3) 



It is well known that, when A > — 1, the global ground- 
state of the XXZ model on the cubic lattice is nondc- 
generate and is a spin singlet |2^, It implies that 

(S?) = in Eq. (J2J). Therefore, entanglement of the two 
spins and the rest part of the system can be character- 
ized by the von Neumann entropy of the reduced density 
matrix (J2J). We have 

E v = —u + log 2 u + — u~ log 2 u~ 

-A + log 2 A+-A-log 2 A-, (4) 

where A^ = W\ ± z. Obviously, the local entanglement 
combines all three correlation functions (SfS?), (SfSif) 
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FIG. 1: The local entanglement (left) of one-dimensional 
XXZ model and its first derivative (right) as a function of the 
anisotropic term A. 



and (S?S?). They behave differently as symmetry of the 
system changes. In the following, we study how the local 
entanglement is affected by them. 

We start with the one-dimensional XXZ model, which 
can be solved exactly by the quantum inverse scattering 
method 31, 32, 331. Its Bethe ansatz solution reads 321 



/ sinli7(Aj 
\sinh7(A :) 
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where the parameter 7 arises from the anisotropic scale 
A, i.e., A = cos 27, and Xj (j = 1, • • ■ , M) are spin rapidi- 
ties, which describe the kinetic behavior of a state with M 
down spins. The regime — 1 < A < 1 is characterized by 
real positive 7 while the regime 1 < A by pure imaginary 
7 with positive imaginary part. From the Bethe ansatz 
equation, the energy spectra and the thermodynamics of 
the system can be obtained from the solution of M spin 
rapidities^- In particular, we find the ground-state en- 
ergy as a function of A. Consequently, the correlation 
functions can be calculated by the Hellman-Feynman the- 
orem. It yields 



(SfS? +1 ) = 



1 dE(A) 
N dA 



(6) 
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For one-dimensional spin-1/2 XXZ model, two critical 
points separate the whole region of A into three different 
phases. These are the ferromagnetic phase for A < —1, 
the metallic phase for — 1 < A < 1 and the antiferro- 
magnetic insulating phase for A > 1. In the ferromag- 
netic region, the ground state is doubly degenerate and 
the local entanglement is zero. It means that there is 
no quantum correlation between two parts of the sys- 
tem. Therefore, we shall concentrate on the region of 
A > — 1. In the Ising limit A — > 00, the ground state is 
approximately a superposition of two Neel states. There- 
fore, the quantum fluctuation is zero, i.e., (SfS? +1 ) = 0. 
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FIG. 2: The concurrence (solid line) and the entanglement 
of formation (dashed line) between two nearest neighbors in 
one-dimensional XXZ model as a function of the anisotropic 
term A. 



Moreover, the longitudinal antifcrromagnetic correlation 
(S* S[ +1 ) = —1/4. Consequently, the local entangle- 
ment tends to unit, i.e. E V (A — > oo) = 1, in this 
limit. It is due to the fact that, as A decreases, the 
effect of hoping process become more and more impor- 
tant. Thus the fluctuation term (S?Sf +1 ) = is un- 
drawn from zero, while the ordering term is suppressed. 
Therefore, the quantum correlation between the local 
part and the rest part of the system will be enhanced. 
On the other hand, when A = 0, the Hamiltonian can be 
transformed into a free spinless fermion model through 
the Jordan- Wigner transformation. In this case, the 
nearest-neighbor correlation functions of the system are 
(SfS? +1 ) = (S?S? +1 ) = -1/27T and (S?S[ +1 ) = -l/n\ 
Therefore, we have E V (A = 0) = 1.3675. Obviously, the 
fluctuation term is dominant, as we can see from the two 
eigenvalue of the reduced density matrix: A + = 0.669 
and A~ = 0.033. By Eq. JI), it produces a relatively 
smaller entropy. In this situation, both ferromagnetic 
and antiferromagnetic couplings are expected to weaken 
the dominant A + . Therefore, when the antiferromagnetic 
coupling is turned on, the competition between ordering 
and fluctuation will lead to a maximum local entangle- 
ment at a certain point. As we can see from Fig. ^ it is 
just the critical point A = 1 where the quantum phase 
transition undergoes. By the Bethe-ansatz result[33|], the 
correlation functions at A = 1 take on values (S?S? +1 ) — 

{S?S? +1 ) = (SfSf +1 ) = (l/4-ln2)/3. Then, the local 
entanglement becomes E V (A = 1) = 1.3759. Moreover, 
the local entanglement is a continuous function around 
the transition point. As A becomes negative, the anti- 
ferromagnetic order is obviously suppressed and the ele- 
ments in pij tends to 1/4. This will also result in a 
larger value of E v . Another interesting result is the sin- 
gular behavior of E v at A = — f , as shown in Fig. [I] It is 
due to the infinite-degeneracy at the critical point. Thus, 
in one-dimensional XXZ model, the quantum phase tran- 
sitions can be well described by the local entanglement. 
However, these properties cannot be identified by the 



ground-state energy, which is usually used in analyzing 
the quantum phase transitions. For instance, at A = 1, 
the ground-state energy is actually a continuous function 
and shows no sign of singularities |4jj ■ 

Clearly, the concept of the local entanglements is quite 
different from the one of the pairwise entanglement, usu- 
ally measured by the concurrence |28j . They describe the 
quantum correlation between different part of the sys- 
tem. However, they are closely to each other. To make 
this point more clear, we show both the concurrence and 
the local entanglement of the present model as a func- 
tion of A in Fig. [21 In the region of < A < oo, as 
we explained above, the competition between ordering 
and fluctuation leads to a maximal value for both the 
local entanglement and the concurrence. As A — > — 1 + , 
however, the concurrence is dropped down to zero and 
the local entanglement is pulled up dramatically. This 
fact is consistent with the monogamy properties of the 
quantum entanglement [34j . On the other hand, the con- 
currence is only valid for the description of entanglement 
between two qubits. Therefore, for fermionic models and 
spin models with spin larger than 1, it is not a well de- 
fined measurement for pairwise entanglement, although 
the negativity |35j may partially characterize it to a cer- 
tain extent. In order to study the role of quantum corre- 
lation played in the quantum phase transition, the local 
entanglement, as the simplest and computable measure- 
ment of the block-block entanglement is a good replace- 
ment, as we showed for the extended Hubbard model 1241. 



III. TWO- AND THREE-DIMENSIONAL 
SPIN-1/2 XXZ MODEL 

Next, we extend our study to two-dimensional anti- 
ferromagnetic spin- 1/2 XXZ model. Since there is no 
exact solution for higher dimensional XXZ model, some 
approximate approaches such as the spin- wave theory [36l 
l37j or numerical calculations on finite lattice (3^, have to 
be applied. In fact, with the help of scaling analysis, the 
results of the exact diagonalization are fully consistent 
with those obtained by the quantum Monte Carlo 39] and 
other analytical method [3rj l37|. By the exact diagonal- 
ization technique, we calculate the local entanglement of 
the two-dimensional model as a function of A on both 
4x4 and 6x6 square lattices with periodic boundary 
condition and show the results in Fig. EI F° r A > f , the 
local entanglement is a decreasing function of A. It is 
due to the fact that enhancement of the antiferromag- 
netic order suppresses the quantum correlation between 
a local part and the rest part of the system. On the other 
hand, for < A < f , it is an increasing function of A, 
as we explained above for the one-dimensional case. Fur- 
thermore, the local entanglement reaches its maximum 
at the isotropic point A = 1. However, unlike its one- 
dimensional counterpart, the local entanglement of the 
two-dimensional XXZ model shows a cusp-like behavior 
around the critical point. Such singular behavior implies 
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FIG. 3: The local entanglement of two-dimensional XXZ 
model as a function of the anisotropic term A for 4x4 (circle 
line) 
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FIG. 4: The representation of E v as a function of A (left) and 
correlation length [ill E^ f^ is in unit of lattice constant) for 
one- dimensional XXZ model in the insulating phase (A > 1) 
with different size. 
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FIG. 5: The local entanglement of one-dimensional XXZ 
model as a function of the anisotropic term A for various 
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FIG. 6: LEFT: The first derivative of the local entanglement 
of one-dimensional XXZ model as a function of the anisotropic 
term A for various system size L — 8, 10, 12; RIGHT: The 
scale analysis of the minimum point of its first derivative. 



the existence of long-range-order, which is absent in the 
one-dimensional case. It is consistent with our previous 
result for the two-dimensional XXZ model 5] . 

As is well known, the ground state of antiferromagnetic 
XXZ model is nondegenerate. One of our main results 
on the one-dimensional XXZ model is that the local en- 
tanglement is a smooth continuous function of A around 
the critical point. In the meantime, the correlation func- 
tions decays by power-law in the antiferromagnetic re- 
gion. Therefore, the local quantities of the system, such 
as energy density and the nearest-neighbor correlation 
functions, are not effected by the those spins, which are 
far away from the local pair. As a result, the finite degree 
of freedom in small system guarantees the analyticity of 
these quantities jijj. In particular, the concurrence is 
continuous 0. For the same reason, the local entangle- 
ment should be also continuous. Furthermore, in Fig. 0] 
we show the scaling behavior of the local entanglement 
for A > 1. We find that the local entanglement does 
not obey a length scaling law. It is rather different from 
behavior of other quantities, such as the spin stiffness in 
the XXZ model [42j. However, in two-dimensional XXZ 



model, the maximum of local entanglement at the criti- 
cal point is sharpened, as the size of the system increases 
(See Fig. |3J). This is due to the existence of long-range 
correlation, which leads to a strong dependence of lo- 
cal entanglement on the system size. In other words, it 
is the infinitely large degree of freedoms in the infinite 
system which introduce the singular behavior of the lo- 
cal entanglement. These results are also consistent with 
the previous conclusions on the concurrence of the two- 
dimensional XXZ model |3, |3(| |33 ■ 

IV. ONE DIMENSIONAL SPIN-1 MODEL 

In this section, we extend our study to a one- 
dimensional spin-1 model, whose Hamiltonian is defined 

as 

H = 53 (sfSf + SfSf + AS?S? - P (Si • 5j) 2 ) , (7) 

where (3 is a real parameter. Differing from the case of 
spin-1/2, the ground-state of spin-1 chain contains fasci- 
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FIG. 7: The local entanglement of one-dimensional bilinear- 
biquadratic model as a function of the 6 for system size L = 6. 



nating surprise which originates from the Haldane gap. 
For the isotropic antiferromagnetic Heisenberg chain with 
integer spin, the existence of this gap in the excitation 
spectrum of the system was first conjectured by Haldane 
|43j |. Then, it was confirmed by the quantum Monte 
Carlo simulation^! and scaling analysis from the exact 
diagonalization 45j . However, as A increases, this gap 
will vanishes around A ~ 1.18. 

To study the role of the local entanglement in connec- 
tion to the critical behavior of this system, we use the 
exact diagonalization method to compute the two-site 
entanglement in its ground state. The results are shown 
in Fig. Obviously, the behavior of the local entan- 
glement is quite similar to the one of spin- 1/2 model. In 
particular, the local entanglement has a maximum at the 
isotropic point, although it is not a quantum transition 
point. A careful scrutiny reveals that the transition oc- 
curring in spin-1 chain is of the Kosterlitz-Thouless type, 
which is quite different from the one for the Heisenberg 
model of spin-1/2. In the latter case, the spin excitation 
spectrum of the system is gapless on both sides of the 
transition point. However, for the system with spin-1, 
the spectrum becomes gapful on one side of the transition 
point. This situation is very similar to the one observed 
for the transverse field Ising model|27]]. Since the first 
derivative of the concurrence is singular and obeys scaling 
law around the transition point of this modelQ, we spec- 
ulate that the same behavior may be seen for the local 
entanglement around the transition point for the spin-1 
Heisenberg chain. Indeed, by taking the first derivative of 
the local entanglement of Hamiltonian Q with respect to 
A, we find a minimum point around A w 1.3 in Fig. 
Moreover, its value varies as the size of the system in- 
creases. By fitting data with respect to the sample size, 
we see that it tends to A ~ 1.18 in the thermodynamic 
limit. 

To explore the effect of parameter /?, we set A = 1 and 



rewrite Hamiltonian J7J as 

H = J2 (cos05i • Sj + sin 9 (§i • S,f) . (8) 

By introducing trigonometric functions, we are also able 
to take the effect of coupling sign into our consideration. 
This Hamiltonian has a rich phase diagram at zero tem- 
perature. It consists of the Haldane phase, trimerized 
phase, and dimerized phase |4q. We now study the role 
of local entanglement around critical points. 

Our results are shown in Fig. \7\ From the figure, 
we find that the local entanglement reaches a local min- 
imum at 9 = 7r/4, which separates the Haldane phase 
and trimerized phase. At 7r/2 < 9 < 57r/4, the ground 
state is ferromagnetic and degenerate. In this case, the 
local entanglement cannot be well defined because the 
thermal ground state comprises all states of lowest en- 
ergy with equal weight. Indeed, a fully polarized state is 
separable. Therefore, if the ground state is spin polar- 
ized one, its entanglement is zero. As a result, sudden 
changes in the local entanglement occur at both 9 = tt/2 
and 9 = 57r/4. They are caused by the ground-state 
level-crossing. However, we do not find any discernible 
structure around the critical point 9 = 77r/4, which sepa- 
rates the dimerized phase from the Haldane phase. Since 
we consider only the sample with L = 6n sites, whose 
ground state has both of trimerized and dimerized order- 
ings, the rapid increasing size imposes further limit on 
scaling analysis. On the other hand, we expect that an 
extremum will appear in its derivatives around 9 = 77r/4 
when the system size becomes large, just like the case of 
A ~ 1.18 in XXZ model. We also find a local minimum 
at 9 = 3tt/2. Around this point, cos# = and the Hamil- 
tonian is reduced to H = — X}(y) (§i ' ^j) • Its ground 
state is nondegenerate, while its first excited state is 3- 
fold degenerate on one side of the transition point and 
5-fold on the other side, the degeneracy is exactly 8-fold 
at the transition point. Moreover, in the previous works 
|46| . it has been shown that the dimerized phase in the 
region 5tt/A < 9 < 3n/2 is gapless. Therefore, a level- 
crossing between the lowest excited states must occur, as 
the one observed for the XXZ model of spin-1/2 at the 
transition point A = 1. Since the minimum in local en- 
tanglement is intrinsically related to the symmetry of the 
Hamiltonian, we believe that the point 9 — 37r/2 is also a 
critical point at which two ordered phases are separated. 

V. SUMMARY AND ACKNOWLEDGMENT 

In the present paper, we study the global phase dia- 
gram of the quantum spin models with either spin-1/2 
or 1 by investigating the local entanglement. We show 
that, indeed, many global properties of the system can 
be derived from such a local measurement. In fact, one 
has observed a long time ago that the original three- 
dimensional image can be recovered from a small piece of 
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holograph, though its resolution is reduced. It is caused 
by the classical interference. Similarly, our findings can 
be understood on the basis of the quantum superposition 
principle. 

We see that, for one-dimensional spin- 1/2 XXZ model, 
the local entanglement shows singular behavior around 
one critical point A = — 1 and takes on its maximum at 
another critical point A = 1. For the two-dimensional 
spin- 1/2 XXZ model, we find that the maximal point of 
local entanglement around A — 1 is sharpened. It is due 
to the existence of long-range order, which is absent in 
one-dimensional case. For the spin-1 XXZ system, the 
local entanglement also has maximum at the isotropic 



point. Moreover, the scaling analysis manifests that its 
first derivative will tend to the critical point A « 1.18 as 
the system size becomes infinite. The rich phase diagram 
for the bilinear-biquadratic model can be almost mapped 
out from the behavior of the local entanglement. Further- 
more, from the singular behavior the local entanglement, 
we find that the point 8 = Stt/2 may be also a transition 
point of this model. This issue deserves definitely further 
investigation. 
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No. 401703) and the Chinese National Science Founda- 
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